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An equation for the quasi-static soliton ansatz depending on an arbitrary 
set of collective variables is covariantly derived on the basis of the variational 
CN \ approach to the method of collective variables. The field configuration and the 

static stresses of the deformed $ 4 - kink that are produced by the excitation 
fNl ' of the internal soliton mode are exactly determined. The kink interaction 

£^ ■ potential at large distances is considered for the example of the nonlinear 

Klein-Gordon system. A general approach to the problem of exactly deter- 
mining the intersoliton potential for the entire set of physically admissible 
CO \ two-soliton configurations is discussed. 
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The dynamic questions in solving some soliton problems are relegated to the back- 
ground, at least, at the initial stage, and the problem of determining the static 
stresses in the system becomes the key point in the investigation. This problem 
directly reduces to finding quasi-static soliton configurations described using a set 
of non-Goldstone collective coordinates. In this context, the problem of determining 
the intersoliton potential in various models can be mentioned. The non-Goldstone 
collective coordinates in these problems are the separation parameter and also the 
^ ■ variables describing the relative orientation of the soliton pair in the internal space 

in multidimensional cases. Furthermore, the non-Goldstone coordinates, which in 
essence are nonlinear amplitudes, describe the dynamics of internal soliton modes. 
The problem of determining the quasi-static configurations in these problems has 
been repeatedly discussed, and the ways to solve it were related to various model 
constructions. In this connection, we mention the investigations of the skyrmion 
interaction begun as early as 1962 ffl, El H, K H. Moreover, the soliton interaction 
was also studied for some other models || |7|, |J, where some important properties 
of intersoliton potentials were established. At, the same time, it is clear that an 
approach based on model semiphenomenological notions does not ensure a compre- 
hensive investigation of the problem, for which a consistent theoretical scheme is 
required. 

Touching upon this problem, we indicate some difficulties related to introducing 
non-Goldstone collective coordinates. For example, the notion of the distance be- 
tween extended objects undergoing mutual deformation acquires its unambiguous 
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visual meaning only if these objects are at sufficiently large distances from each other. 
Therefore, a consistent theory must not only describe the process of deformation of 
solitons as they approach each other but also answer the question of interpreting the 
ambiguity in determining the intersoliton distance at the stage where the solitons 
lose their identity as a result of this deformation. The difficulties in constructing a 
theoretical scheme suitable for describing intersoliton interaction that are due to the 
above-mentioned factor have been extensively mentioned in the literature and were 



most distinctly stated in ||]. Some simple considerations were later presented in [10 
in favor of the universal variational approach to the method of collective variables, 
which permits stating various problems of soliton dynamics in the framework of the 
Hamiltonian procedure. In ||10|| , the general equation for quasi-static configurations 
determined by an arbitrary set of collective coordinates was also derived. It is es- 
sential that the resulting equation is invariant under reparameterizing the collective 
variables, which ensures that the physical quantities are independent of the way in 
which the intersoliton distance is defined. 

In Sec. |2|, we present a new covariant method for deriving an equation for 
the quasi-static ansatz using basic notions related to Riemannian manifolds. The 
following two sections are devoted to elaborating solution methods for this equation 
in some specific cases. In Sec. |3], the problem of determining the configurations and 
static stresses of the <3> 4 -kink whose deformation is produced by the excitation of the 
internal soliton mode is considered. This problem can be solved exactly. Similar 
problems were considered in many papers. The corresponding references can be 



found in [11, 12 1 . It is unlikely that there is a model with a nontrivial potential that 



admits an exact solution of the problem of determining the intersoliton interaction. 
A solution of this problem that allows calculating an asymptotic expression for the 
intersoliton interaction can most likely be found for an arbitrary model. An example 
of a calculation of this kind is given in Sec. f|, where the equation for the quasi- 
static ansatz is used to determine the kink interaction for the nonlinear Klein-Gordon 
equation in the domain of large intersoliton distances. We stress that, apart from 
independent interest, knowing the intersoliton potential in the asymptotic domain 
permits introducing a physically acceptable gauge in the numerical analysis of the 
equation throughout the range of intersoliton distances. In Sec. [5], we discuss the 
investigation results and indicate the main elements of the suggested theoretical 
scheme. 



2 Equation for the quasi-static ansatz 

We consider the question of finding the quasi-static ansatz for a soliton system 
described by an arbitrary set of collective variables. Let 

be a nondegenerate Lagrangian quadratic with respect to velocities that admits soli- 
ton solutions, let q l be the field coordinates, and let i = {i,x} be the set of discrete 
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indices and spatial coordinates. It is convenient to interpret the basic configuration 
space as a Riemannian manifold with a metric generated by the kinetic term. In 
this case, if Q are collective variables, then the quasi-static ansatz q\{Q) describing 
the deformed soliton state can be interpreted as a submanifold parameterized by the 
variables Q. We assume that this ansatz is determined by the stationarity condition 
for the Hamilton function defined on q l c (Q) for the system in question relative to 
arbitrary infinitesimal variations Sq l c (Q) orthogonal to the manifold surface, i.e., 

«* = 0, (1) 

where the variations 5q l c (Q) satisfy the conditions 

9M^Sq k c (Q) = 0. 

To cancel the variations in Eq. ([[]), we introduce the operator V of projection onto 
the subspace tangent to the submanifold q l c {Q) at the point Q 

n^wg^g, v^v, ( 2 , 

where 



dql ,<9g, 
:9ik{q c 



is the metric on the submanifold q l c {Q) and g®® is the matrix inverse to (?qq. Because 
the variations Sq l c (Q) in the orthogonal directions to q l c (Q) are arbitrary, Eq. ([[J) 
can then be brought to the form 

il-V) — = 0. 

oq c 

Substituting expression (|2]) for the operator V in this formula, we obtain 

dH(q c ,0) _ dq k c QQ dE(Q) 

dqi -9Mqq9 -qq-, W 

where E(Q) = H(q c (Q),0) is the energy of the quasi-static configuration. By the 
definition of the tensor g®®, the right-hand side of Eq. (|3]), which takes the exis- 
tence of static stresses in the collective subsystem into account, is invariant under 
an arbitrary nondegenerate change of variables Q — > Q'(Q). This indicates that 
the variational procedure under consideration is covariant. In this connection, we 
note that in relation to the problem of determining two-soliton configurations, for 
example, Eq. exactly describes these configurations in the domain of distances 
where the deformed solitons lose their identity. The choice of the gauge only deter- 
mines the specific form of the parameterization describing the configurations. For 
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convenience, it can be required that the separation parameter coincide with the nat- 
ural intersoliton distance if this distance notably exceeds the soliton size. We note 
that if the collective coordinates are Goldstone variables, i.e., if the coordinates Q 
coincide with the degeneration parameters of vacuum solutions, then the energy E 
does not depend on Q, and Eq. (H) becomes the equation dH(q c ,0)/dq l c = for 
static solutions. 

3 Static stresses of the deformed <£> 4 -kink. 

We use Eq. ([|) to determine the configuration of the deformed kink in the $ 4 model. 
The Lagrangian of the model and the corresponding static kink solution u c (x — X) 
defined to within the spatial coordinate X have the forms 



The eigenfunctions of the operator C describing linear fluctuations on the back- 
ground of the kink solution include the mode cr(x) belonging to the discrete-spectrum 
states, 



£a(x) = 3a(x), a(x) = — , £ = -— + 4 - 6 cosh" 2 x. (5) 



In describing nonstationary states, it is usually assumed that the $ 4 -kink possesses 
an internal degree of freedom whose small oscillations can be identified with the 
mode a(x). Therefore, apart from the translational coordinate X, the set of collec- 
tive coordinates describing the kink state must also include the variable r related 
to the excitation of the internal mode of the $ 4 -kink. For this set of collective coor- 
dinates, we can expect that the directing vectors <&' c = d§ c /dx and $ CT = <9$ c /<9r 
constructed on the solutions of Eq. (fj) are orthogonal in the tangent space, 



Also, taking into account that the energy is independent of the cyclic variable, we 









(7) 



The right-hand side of Eq. (|7|) is invariant under the reparameterization r 
If the collective variable is fixed using the gauge 




- < $ 2 T >^ dE/dr = 1, 
then Eq. (|7]) becomes the nonlinear diffusion equation 





(9) 
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Equation (g) was investigated in Hl3fl , where its exact solutions were found. In 
particular, the solutions of the type of a deformed kink that tend to +1 as x — > +00 
can be written in the form 



e 2x 



^■ r ) > + l ± ^ (10) 

Integrating, we can easily verify that these solutions satisfy orthogonality condition 
(Bp and gauge condition (|8|). We now note that the physical domain of small de- 
formations, where solutions (|10D only slightly differ from static kink solution (f|), 
corresponds to r — ► 00. In other words, the collective variable r defined by gauge 
condition (|j) is inconvenient from the physical standpoint. It plays only an inter- 
mediate role as it permits representing Eq. in the form of a nonlinear diffusion 
equation. It is therefore advisable to combine the two solutions in ( |I~0D to form one 
solution, 

$c(X ' g) = e 2x + x _ ^ x ~ Uc ^ + qCT ^ + ( U ) 

where the new collective variable q, which takes both positive and negative values, 
has the meaning of the amplitude of the internal mode under small deformations 



1C|. The configuration energy E(q) = if($ c ,0) determining the static stresses of 



the deformed kink has the form 



E{q) = t + ih? ~ (1 _V) 3/2 arccos ^ ^t + l q2+ 0(g3) ' 1 = 

Here, the first term is the kink mass, and the quadratic term with respect to q 
describes the potential of a harmonic oscillator with the internal mode frequency 
uj 2 = 3. The range of the collective variable where the energy remains real is confined 
to the interval —a/2/3 < q < + a/2/3. Solution (p|) has no singularities with respect 
to the spatial coordinate in this range of q. We also stress that the choice of the 
physical collective variable tending to the amplitude of the internal mode in the 
limit of small deformations is ambiguous. Namely, it is determined for convenience 
and is connected with explicit form ( PH| ) of the solution for the deformed kink in the 
situation in question. 



4 Kink interaction. 

An important role in studying the properties of soliton states is played by the so- 
lutions of Eq. (£jp that are defined by the asymptotics & c (x) — » — 1 as x — > ±00. 
These solutions describe the kink-antikink configurations and permit establishing 
the properties of the intersoliton interaction. In this case, the collective variable 
relates to the separation parameter r, 2r is the intersoliton distance. The excita- 
tions of internal modes are not taken into account. Unfortunately, it seems that 
exact solutions of Eq. (|j) with given boundary conditions have not yet been found 
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Hfl . Approximate solutions in the domain of large distances can be obtained in the 



general form if the corresponding equation 

+ U'{* e {x)) = ^ < $L >" 1 ^ (12) 

for the nonlinear Klein-Gordon model is considered instead of ((7|). Here, U($) is 
the potential of the model. If the origin is placed at the center of the bisoliton 
pair and its symmetry $ c (x, r) = ±$ c (— x, r) is taken into account, then the entire 
configuration can be described by considering the domain x ^ 0. The even functions 
describe the kink-antikink states and the odd ones describe the kink-kink states (in 
the case of a periodic potential U (<&)). If the intersoliton distance is large, then to 
the left of the point x — 0, the ansatz <& c (x, r) only slightly differs from the one-kink 
solution u c (x + r) of the static Klein-Gordon equation 

-u" c (x) + U'(u c (x)) = 0. 

Therefore, representing the ansatz <& c (x, r) in the form 

§ c (x, r) = u c (x + r) + i](x + r, r), x ^ 0, (13) 

where r](x + r, r) tends to zero as r — > oo, and taking into consideration that in the 
domain of large values of r, 

oo 

< $ 2 cr >^2 dxu'?{x) = 2M, (14) 



where M is the kink mass, we derive an equation for r)(x) from fll2"D , 

1 r)E rfi 
C7]{X) = 2MW <(X) > * <r ' £ = -q^ + U"Mx)). (15) 

Here, dE/dr is a small parameter in the domain of large values of r. Multiplying 
Eq. ( |i~5"D by the zero mode u' c (x) of the operator £ and integrating by parts from 
— oo to x , we can reduce the order of the equation. As a result, we obtain the linear 
first-order equation for 77(2;) 

X 

-u' c (x)r,'(x)+u'>(x)ri(x) = M(x) = J dx'u'^x'). (16) 

—00 

The solution of Eq. (|IED can be written as 

1 V(r) 1 dE f , ,M(x') . 
^) = <(*) 373 + 77773 / ^777777 ■ ( 17 ) 



<(r) 2Mdr J u'^(x') 
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The function i](r) playing the role of the integration constant in (fL7j) is found from 
the boundary conditions. For the kink-antikink and kink-kink pairs, these conditions 
have the respective forms $' c (0, r) = and $ c (0, r) = 0. Using these conditions and 
Eq. (|16D, we obtain the following expressions for the boundary values r]{r) and rj'(r) 
in these two cases: 

1 (M{r)dE , 2/ 



V (r) = -u c (r), v{r) = [^fj^ ~ u c (0 ) ( ' > 



for the kink-antikink pair and 



r,(r) = -u c (r), rf(r) = + <VK(r) ) (19) 



for the kink-kink pair. Solution fll7|) describes the distortion of the kink under 
the effect of the second soliton when the distance between the kinks substantially 
exceeds their sizes and depends on the small parameter dE/dr, which determines 
the interaction of the solitons. We note that formulas ( JTT| - |T§D include terms of higher 
order than the leading asymptotic terms with respect to r. However, because the 
function i](x) contains competing spatial exponentials, it is inadvisable to select 



these terms in considering the field variable rj(x). Formulas ([17l - |T9|) should therefore 
be used to correctly determine the leading approximation term for the integrated 
characteristic, the interkink potential E{r). Because {l/2)u'^ = U(u c ), substituting 
(|l3l) in the formula for the energy E(r) = H($ c , 0), expanding with respect to rj(x), 
and integrating by parts with regard to Eq. (|T6|) result in 

oo r 



/l 9E f 
dxu?{x) + —-^ J dx V (x)u' c (x) 



+ r](r)(2u' c (r)+v'(r)). (20) 
Substituting expression ( |P7D for rj(x) in ( p0|) and using boundary values fll8|) and 



(O) for the kink-antikink and kink-kink pairs, we see that the linear terms with 



respect to dE/dr cancel. In this case, relation ( P0|) becomes 



1 dE^ 2 



E(r)-2M = A(r) + B(r)(— — ) , (21) 



where the functions A(r) and B(r) are given by the formulas 

3 °° 

A(r) = -^}f\-2 I dxu'^(x), 



M 2 (r) f , M 2 {x) 
B(r) = / \ / \ + dx- 



u' c (r)u"(r) J u'c( x ) 
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for the kink-antikink pair and 



A(r) 



M c( r )~frr + 2 / dxu c (x)u"(x) , 
u' c {r) > 



B(r) 



dx 



M 2 (x) 
u' 2 {x) 



for the kink-kink pair. The asymptotic expression for the kink solution has the form 

u c (x)\ x ^ ±00 ^u ± =F aexp(=Fmx), 

where m = \/U"(u±) is the pion mass. In a case of an antisymmetric kink-kink 
configuration the vacuum constant u + is equal to zero. Taking into consideration 
that M(oo) = M and M(— oo) = and calculating the leading approximation terms 
with respect to e~ mr for the coefficients A(r) and B(r), we see that the functions A(r) 
vanish and B(r) coincide to within the sign. Therefore the leading approximation 
for Eq. fl2ip has the form 

E(r) -2M = T^r~ 2 e 2mr ( ^ 
om^cr \ or 

whence we find the interaction energy of the bisoliton pair at large distances, 

E(R) = 2M =f 2maV 2mr . (22) 

This formula describes attraction (repulsion) for the kink-antikink (kink-kink) con- 
figuration. The expression for the energy in ( |22] ) in the case of the $ 4 and sine- 
Gordon systems was found in || using an approximate equation of type ([12]) whose 
right-hand side was imitated with a phenomenological delta-shaped source [I5|. 



To determine the soliton configuration and interaction throughout the range of 
distances, the gauge should be fixed, and Eq. (|TJ) should be solved numerically. 
We note that the representation of the solution in form ( |T3D partly fixes the gauge 
in the domain of large intersoliton distances. Formulas ([14]) and (|22]) obtained using 
expansion ([13]) can be used to introduce the natural gauge 

l9 i BE 2m 2 a 2 9mr . , , 

< < > ' = (23) 

ensuring the coincidence of the doubled parameter r with the true intersoliton dis- 
tance in the asymptotic domain r — ► oo. After gauge (|3|) is fixed, Eq. ( |i~2|) takes 
the form of a nonlinear diffusion equation, 

-*» W + ^x)) = ± ^ e --^M. (24) 
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Equation (^4|) can be used to find the corrections to potential (22) in the asymptotic 
domain for the given gauge, and its numerical solution permits reproducing the entire 
set of physically admissible two-soliton configurations. We stress that solving this 
equation in the case of the kink-antikink configuration corresponds to describing 
diffusion in the direction of positive "time" — r. This process is stable, and the 
mutual approach in the kink-antikink pair terminates with its "annihilation" under 
which $ c (x) goes into the vacuum constant f!4 |. Actually, the calculation of the 



evolution of the kink-antikink pair performed in 0] in relation to the investigation 
of the nonlinear diffusion equation is an example of rigorous determination of the 
entire set of two-soliton configurations in problems related to studying intersoliton 
interaction on the basis of Eq. ([|). 

Constructing the kink-kink configuration reduces to studying diffusion in the 
direction of negative "time" r. In the domain of large intersoliton distances, the 
presence of nonlinearity in Eq. (|24|) stabilizes the instability characteristic of this 
procedure in the case of the linear theory. The stabilization weakens during the 
mutual approach, and the limiting field configuration is determined by the point 
where the stability of the process terminates. 



5 Conclusion 

We have covariantly derived an equation for a quasi-static ansatz on the basis of 
the variational approach. A characteristic feature of this approach is the invariance 
under reparameterization with a change of collective variables. The cyclic variables 
or the variables related to the degeneration of vacuum states are introduced in the 
standard way. The arbitrariness in the choice of the non-Goldstone variables can 
be eliminated based on some physical considerations in the domain where these 
variables can be given an immediate physical interpretation. Otherwise they are 
chosen for convenience. 

As an example, we have considered the static deformations of the $ 4 -kink re- 
lated to excitation of the internal mode. This problem admits an exact solution. 
The configurations of the deformed kink, its energy, and the physically admissible 
range of the amplitude of nonlinear oscillations have been determined. Determin- 
ing the intersoliton potential for the nonlinear Klein- Gordon equation has also been 
touched upon. Linearizing the equation for the quasi-static ansatz in the domain 
of large distances permits establishing the deformation of the kinks and the asymp- 
totic character of their interaction, which allows introducing the natural gauge and 
representing the equation for the quasi-static ansatz as a nonlinear diffusion equa- 
tion suitable for numerical determination of the entire set of physically admissible 
two-soliton configurations 

In relation to determining the intersoliton interaction, we stress that although 
two-soliton configurations have been comprehensively studied using various qualita- 
tive approaches for the majority of actual models, investigations on the basis Eq. 
([H) permit including the exact effect of soliton deformation at small distances, i.e., 
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in the domain most sensitive to the defects characteristic of artificial field construc- 
tions. This is particularly important, for instance, in the Skyrme model because 
the nucleon-nucleon potential in this domain has an insignificant minimum forming 
the nuclear coupling energy. Although Eq.(|3]) does not take the spin effects into 
account, the study of skyrmion-skyrmion interaction on the basis of this equation is 
a useful stage in the procedure of determining the nucleon-nucleon potential. 

The author is grateful to P. 0. Mchedlov-Petrosyan and V.V. Gann for a number 
helpful discussions. 
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